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l/y ' Abstract 



We consider the relativistic generalization of the harmonic oscillator problem by addressing dif- 
ferent questions regarding its classical aspects. We treat the problem using the formalism of Hamil- 
tonian mechanics. A Lie algebraic technique is used to solve the associated Liouville equations, 
yielding the phase space evolution of an ensemble of relativistic particles, subject to a "harmonic" 
potential. The non-harmonic distortion of the spatial and momentum distributions due to the in- 
trinsic non-linear nature of the relativistic contributions are discussed. We analyze the relativistic 
dynamics induced by two types of Hamiltonian, which can ascribed to those of harmonic oscilla- 
tors type. Finally, we briefly discuss the quantum aspects of the problem by considering possible 
strategies for the solution of the associated Salpeter equation. 
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I. INTRODUCTION 



The relativistic harmonic oscillator (RHO) is a topic of fundamental importance in 
Physics. Although widely discussed in the past (see, for example Ref. it still contains 
elements of interest, deserving careful comments, both in classical and quantum mechanics. 
The spinless Salpeter equation with a quadratic potential cannot be solved in analytical form 
[2! and the corrections induced by the relativistic effect in the phase space evolution of an 
ensemble of particles ruled by a harmonic potential have not been discussed yet. This point 
is raised by a genuine practical necessity, emerging in charged beam transport in magnetic 
systems. In particular, in the case of propagation through quadrupole lenses, the beam is 
ruled by a potential of harmonic oscillator type jsl. It may therefore happen that, in ab- 
sence of a clear distinction, phase space distortions due to relativistic corrections could be 
misinterpreted and ascribed, for example, to quadrupole aberrations. The dynamics of the 
relativistic harmonic oscillator is, indeed, not harmonic at all. The relativistic nature of the 
problem induces non-linear corrections which determines a genuine anharmonic behavior. 

In this paper we will consider the classical aspects of the RHO and we will treat the 
problem ab initio. Just to give an idea of the problems we are going to treat, we consider 
the case of a relativistic particle subject to a linear potential whose Hamiltonian is 

H = cy^f+ljn^ - Fx, (1.1) 

and the associated equations of motion 

dH pc dH 

a; = — — = — , =, p = —^1— = F = rnoa, (1.2) 

dp ^p^ + imo c)2' dx ^ ' 

are readily integrated as 

x{t) = xo+ (Vi+TP - J I + Ul] -, U{t) = ^ = Ho + — . (1.3) 

\ ^ J a itiqc c 

The non-relativistic limit of these equations is given by 



x(^^\t)=xo + l 



no + — j -n^ 



n(^^)(t) = no = ^ = 7o/3o, (1.4) 
mo c 



{(3 = f/c, 7 = l/yl — /3^) and is recognized as the classical equation for the uniformly 
accelerated motion. 



2 



It has been noted by Harvey {4] that the derivation, in an unambiguously Lorentz invariant 
fashion, of the Hamiltonian of a relativistic particle subject to a scalar potential, hinges on 
the relaxation of the requirement that the rest mass be a constant. The procedure proposed 



m 
in 



Ref. I4 



J] is reminiscent of the Lorentz covariant theory of gravitation (see last reference 
^), where the rest mass is potential dependent. By adhering to this point of view, an 
alternative relativistic Hamiltonian with a linear potential can be written, according to the 
prescription 



H — [pc) = {mc ) , m = rrio — (1.5) 



Fx 

c 

where it is evident that the potential redefines the particle mass, and the following Hamil 
tonian is obtained 



H = v/(pc)2 + (mo c^-Fxy. (1.6) 

The physics described by this Hamiltonian totally differs from that of Hamiltonian (11. ip . 
For example, the latter can be exploited to study the equation of motion of a relativistic 
charged particle subject to an external electric field, while Eq. (II. 6p could be viewed as that 
relevant to a relativistic charged particle in a constant magnetic field^. 

To better appreciate the differences, we make a comparison between the solutions of the 
respective equations of motion, which in the case of the Hamiltonian (II. 6p . read 

. Pc^ . mo - Fx 

a; = -^, P = F (1.7) 

and can be integrated to get 

2 

x(t) = ^ {l + 7o (v^r^sine-5 cose)}, (1-8) 

with 

^ at e 1 ^ ^ 

7o 7o ^ / 

The motion, as also evident from Eq. (II. 7p . is essentially harmonic with a characteristic 

frequency u = 0/(70 c). It is interesting to note that, the solution (including the velocity) 

can be written in the form of a rotation, namely 

sin ^ — cos ^ 
cos e sin e 




)l 









[1.9] 



^ The relativistic Hamiltonian for a charged particle with charge q and mass too moving in a magnetic field 
with intensity B directed along the z-direction can be written as H = c ^Jp^ + {py — qB x)^ + (toq c)^, 
mathematically equivalent to Eq. (|1.6p . 
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where 

(1.10) 

and, thus, Xq = Vl — 5^, = 6. From this result we see that the particle motion is 
essentially oscillatory, and for short times {t ^ c-fo/a) the non-relativistic equations of the 
uniformly accelerated motion is obtained. 

In the forthcoming sections we will discuss analogous problems in the case of the harmonic 
potential. We will address our analysis by studying the Liouville equation describing the 
evolution of phase-space distributions of an ensamble of relativistic particles. 



^(0 
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II. THE LIOUVILLE EQUATION FOR THE RELATIVISTIC HARMONIC OS- 
CILLATOR 

Let us now consider the following Hamiltonian 



H = cy/p"^ + {mocy + V{x) (2.1) 

which rules the relativistic evolution of a particle with mass rriQ subject to a potential 
V{x). The solution of this dynamical problem is accomplished by solving the associated 
Hamiltonian equations, that introducing the notation Z = (x,p)"^, can be written as 

Z = HoZ, (2.2) 

where 

HoF( )- \H F\-^^^^ (2 3) 

dp dx dx dp 

The solution of Eq. (12. 2p can be formally expressed as the following Lie series 

Z = Unit) Zo = y2- {Hoy Zo (2.4) 

n=0 

where Unit) = e*^° is the Hamilton evolution operator, and the following notation has been 
used 

(i/o)OF = F, {Ho)F = {H,F}, iHoYF = {H,{H,F}}, ... (2.5) 

This procedure is reminiscent of analogous methods adopted in quantum mechanics to treat 
the time-dependent Schrodinger equation . The problem is not in general amenable for an- 
aljd;ical solutions. For example, in classical (non-relativistic) mechanics exact solutions can 



be obtained only in case of quadratic Hamiltonians. In the relativistic case even quadratic 
potentials are not amenable for an exact solution. In these situations, the search of solutions 
of approximate nature is the only way to address the problem. The symmetric exponential 
split rule has been proven an efficient tool. According to it, we can write 

^tHo ^ ^it/2)H,o ^tH,o ^m)H,o ^ ^(^3) (2.6) 

where 



Hi = c ^v" + (^0 c)2, = V{x). (2.7) 

The solution of the problem is then obtained by a repeated application of the evolution 
operator calculated for short time interval steps of duration 6t = t/N, namely 



e''''°Zn (n = 0,l,...,iV-l). (2.8) 



Since we are interested in particle ensemble distributions ruled by relativistic Hamiltoni- 
ans, we can adopt a complementary procedure based on the solution of the Liouville equation 
associated with the Hamiltonian H. It describes the evolution of the particles phase space 
density distribution p{x,p;t) and reads (t is time in the laboratory frame) 

d 

— p(a;, p;t) = L p{x, p; t) (2.9) 

whose solution can formally be written as 

p{x,U-t) = UL{t)p{x,U-0) UL{t) = e'^ (2.10) 

where the Liouville operator L is given by 

'dH d dH d 



(2.11) 



dx dp dp dx 

The use of well-known techniques from operator calculus [5| will be exploited in the forth- 
coming section to discuss the problems concerning the relativistic phase space conservation. 
Here we just remind that the function p{x, p; t) transforms as a Lorentz scalar and its integral 
over the phase space variable is normalized to unity. 

In the case of free relativistic particles the Liouville operator is^ 

n ' (2,12) 



VTTTp dx 



^ We have rescaled the phase spaee eoordinates and exploited the adimensional momentum 11 introduced 
in the previous seetion. 
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and, thus, the action of the associated evolution operator on the initial distribution reduces 
to a shift of the spatial coordinate, namely 



p(x, 11; t) = exp < —ct 



n 



d 



p [ X — ct ■ 



VTTTP dx 

n 



n;0 



The Liouville operator associated with the Hamiltonian ( 11. ip is 



(2.13) 



L = -c- 



n 



d ad 



^/TTTP dx c dn 
and the evolution of the corresponding phase space distribution function writes 



(2.14) 



p{x, Il;t) = p \ X + 



n 



at 



„ at 
n- — ;0 



(2.15) 



We will use this equation as a benchmark of an approximate solution based on the split 
operator method, then applied to the case of more complicated potentials for which no 
exact solution is available. We start decomposing the Liouville operator in two parts 

c n d 



- ad 



(2.16) 



Since Li and L2 do not commute, the evolution operator cannot be naively disentangled, 
the use of the same procedure leading to Eq. (12. 8p . yields 



Pn+l — Ul pn 



(n = 0,l,...,iV-l) 



with 



Ul = exp ( ^Lij exp{6t L2) exp (^Li 



(2.17) 



(2.18) 



By repeatedly acting with this operator on the initial function po = pi^o, Hq] 0), and taking 
into account the shift of the coordinates x and 11 induced at each step by the exponential 
operator, we end up with 

p„ = p(x„, n„; 0) (2.19) 

where 



Xn "^n— 1 



n„ 6t 

2c 



1+ u-^6t 



: C6t, 



2c 



n„ = n„_i--5t. (2.20) 

c 



An advantage of the symmetric split is that of providing at any step a norm preserving 
distribution. The comparison between the exact solution and the one obtained with this 
technique has been completely satisfactory, and therefore the method can be safely applied 
to the case of a quadratic potential. 

We consider the relativistic Hamiltonian 

H = c^/p^^ + {mo cy + l-kx^ (2.21) 

which for future convenience will be written as 

H = mQ^ (^l + n2 + ^ rf^ (2.22) 

where 

r7 = -x, — • 2.23 

C V "^0 

By setting \ = VLt, the associated Liouville operator reads 

r/— . (2.24) 



and the use of the symmetric split method yields the iterated solution 

p„ = p(r/„,n„;0), (2.25) 

with {5\ = n 5t) 

r,n = Vn^i - F' {lin-i + 5A, F(n) = Vl + n2, (2.26) 

and 

n. = n„_i + ^" + ^"-^ 5\. (2.27) 

Fig. [1] show the time evolution in phase space of an initially Gaussian distribution. 

A global idea of distortion induced by relativistic effects in the distribution is offered by 
Fig. 121 where we have reported the 3-dimensional plot after 2400 steps. The distortion is 
more significant on the tail of the distribution, where the momenta are larger and, thus, 
the relativistic effects are more effective. The evolution differs significantly from the corre- 
sponding non-relativistic case, where the distribution remains Gaussian and the evolution 
of the phase space contour plots are realized through continuous deformations of equal area 
ellipses. The distribution of spatial and momentum variables, defined by 

/OO /"OO 
dnp,, i?n(n)= / dr//7„, (2.28) 

-OO J — OO 



FIG. 1: Phase space contour plots after n steps 5A = 5 x 10 (a) n = 0; (b) n = 1000; (c) 
n = 1400; (d) n = 2400. 

are shown in Fig. [3l 

Let us now consider an analogous (11. 6p for the Hamiltonian of relativistic harmonic 
oscillator, namely 



We discuss this example in order to test the the efficiency of the method on a more com- 
plicated mathematical structure than the previous case, and to deal with a relativistic 
Hamiltonian containing higher order non-linearities in the spatial coordinate. Such a model 
Hamiltonian could, in principle, be exploited to study the relativistic motion of electrons in 
non- homogeneous magnetic fields j^. The associated Liouville equation reads 






(2.30) 



and the use of the symmetric split method yields the following iterated solution 



Pn = piVn, n„; 0) 
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FIG. 2: Three-dimensional phase space density after 2400 steps. 




FIG. 3: The distributions defined in Eq. (|2:28]) after n steps 6\ = 5 x 10"^: (a) n = 1000; (b) 
n = 1400; (c) n = 2400. 

with 

Vn = Vn-1 - I ^ SX (2.32) 

SW A 1 2 



n„_i + r]n-i Yj + ^1 + 2 ''"-I 

and n„ given by Eq. f l2.27p . In Fig. Hlthe coordinate and momentum distributions for the 
Hamiltonians fl2.2ip and fl2.29p are compared. 
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FIG. 4: The distributions of coordinate and momentum at the step n = 2400 for Hamiltonians 
(|2:2T]1 (curve (a)) and ([2:29]) (curve (b)). 



III. CONCLUDING REMARKS 

Even though we have insisted on the term relativistic "harmonic oscillator" , the relevant 
phenomenology is by no means harmonic. The interplay between energy and mass in rela- 
tivistic dynamics induces non-linear terms that are responsible for anharmonic contributions 
in the equations of motion and can be, in principle, analytically accounted for by elliptic 
functions. This is not the case for the equation obtained from equation of motion associated 
to Hamiltonian f l2.2ip . namely 

(3,1) 



If we limit ourselves to the lowest-order corrections, we can expand this equation to get the 
following Duffing-type oscillator equation 

n" = -n + ^ (3.2) 

ihat is accurate up to 0(/3^), and whose solution is written, in terms of Jacobi functions, as 

3 

n = V2(l-^2)cn (^a A; (3.3) 

for the initial conditions 

A = o, no = n(o) = ^2(1-^2), n'o = ^nu=o=o. (3.4) 



The oscillation period of the function (13. 3 p is given by 

ng \ . , 2 vr 



T = K 



10 



i + 7n^ 1 



2n2 -4 



(3.5) 



where K is the complete elhptic integral of the first kind. On the r. h. s. we have reported 
the lowest order corrections, from which we may argue that even small kinetic energies (i.e., 
few keV for electrons) can induce corrections to the period of the order of 10%. In Fig. Elwe 
have reported a comparison between the solutions of Eq. (13.1 p obtained by an integration 
with conventional means and via the symmetric split method. The results coincide and 
give complete confidence on the reliability of the previous analysis. On the other side 
the evolution operators associated with the Liouville equations or with the Hamiltonian 
counterparts are linked by the relation Unit) = UL{—t) (see Ref. 8| for further comments). 




-1.0 H — • — • — I — \ — • — I — • — \ — • — I — • — \ — I — • — • — \ — I — • — • — I — • — • — • — I 

6 10 15 20 25 30 

A 



FIG. 5: Solution of Eq. (j3.ip (Hq = 0.9, IIq = 0) obtained with standard numerical methods and 
with the Lie series. The results coincide (blue continuous line). The dot-dashed line refers to the 
ordinary harmonic oscillator. 

In this paper we have addressed the problem of phase space evolution for a relativistic 
system. We have used the Liouville equations, without taking any care about the transfor- 
mation properties of the associated distribution. For this particular problem we have two 
different type of transformations: 

i) those associated with the iterative scheme deriving from the sjTiimetric split decom- 
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position of the evolution operator, 

ii) the Lorentz transformation properties of the Liouville distribution. 

As for the first point, it is worth noting that they can be considered canonical transfor- 
mations, for which, at each time step, we redefine the canonical coordinates (?7„, n„) — )■ 
(?7n+i, n„+i). It can be checked by direct computation it is verified the condition 



J 



drjn 


drjn 






dUn 


dlin 







(3.6) 



that ensures the preservation of the norm of the distribution. 

Although the Liouville distribution is a Lorentz invariant (see Refs. 
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JjJ), its spatial 



and momentum distributions do not possess this property. In particular, we can introduce 
the Liouville density current, defined as (77' = dry/dA) 

/oo 1*00 
dnp(r/,n;A), I{V,\)= dnr/'p(r/,n;A) (3.7) 
-00 J —00 

transforming as a four- vector in Minkowski spacetime and satisfying the continuity equation 

dr,I{v,X) + dxS{v,X)=0. (3.8) 

We limited ourselves to a one dimensional oscillator, but the previous conclusions hold also 
in the three dimensional case. We have checked numerically the correctness of the continuity 
equation (13. 8 p and in Fig. [6] we report the spatial and temporal parts of the Liouville density 
current for the harmonic oscillator. 

In this paper we did not consider the problem from the quantum point of view, by taking 
into account the evolution of the Wigner distribution function, which can be considered the 
quantum analogous of the Liouville phase space density function. This problem has been 
considered in the literature (see Ref. 12[ and references therein), with particular emphasis 
on gauge invariant forms. The analysis we have developed here (in particular, the use of 
the split technique method) can also be extended to this case but completely new scenarios 
are open, mainly with reference to the equation ruling the evolution of the Wigner phase 
space distribution. The use of the Wigner- Moyal phase space operator representation [l3| 
determines additional computational problems, which deserves more than a paragraph in 
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FIG. 6: Liouville density current for the relativistic harmonic oscihator n steps (^A = 5 x 10 ^: a) 
n = 1000; b) n = 1400; c) n = 2400. 

a concluding section. It is our intention to consider the treatment of this problem in a 
forthcoming investigation. 

A point that we can address, without taking too much space, is the solution of the rela- 
tivistic Schrodinger equation for a spinless quantum relativistic particle under the influence 
of constant electric field. From the mathematical point of view, the equation we should deal 

n 

with is the following Salpeter equation 



The same equation in the Fourier conjugated space reads 
where 



The solution of Eq. fl3.10p writes 



dee-'''«vl/(^,r). 



\l/(?7, r) = exp <j -i j dx V 1 + (?7 - A x)^ r ^(^7 - ^ 



(3.9) 



(3.10) 



(3.11) 



(3.12) 
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and therefore its spatial counterpart can be written in terms of the following integral trans- 
form 

^(e,r) = -|= j dr/exp|z (v^' dx v^l + (r/ - A x)^^ } ^(r/ - A r). (3.13) 
In the case of a Salpeter equation with a quadratic potential, i.e 



(3.14) 



by applying the symmetric split technique, and working in the conjugated space, we can get 
the solution in terms of the following recursion 



exp \-i I ^l + riA foo c ^ 

vl/„+,(r/, r) = ^ / da exp \ --j-^ + 1 V^V^ \ r) 



2ViTTBT 



(3.15) 



(\I/o('7!t) = ^('7;0)) where to describe the action of the transformed quadratic potential 
on the successive steps of the iteration we have used the following imaginary Weierstrass 
transform 



e'"^'^ f{r]) 



1 



daexp<!-^^^^[. /(a, 



(3.16) 



This approach can be complemented with the generalized transform methods developed in 



Refs. 



14| and 15|. As it will be shown elsewhere, this procedure can be efficiently used to 



study evolution of the wavefunction of a relativistic spinless particle. 
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